The static response and the Local-field factor of the 2-D electron fluid. 
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The static local-field factor (LFF) of the interacting 2-D fluid is calculated nonperturbatively using 
a mapping of the quantum fluid to a classical Coulomb fluid [Phys. Rev. Let., 87, 206 (2002)]. 
The LFF for the paramagnetic fluid differs markedly from expectations from standard perturbation 
theory. Our result for the LFF has a maximum close to 3/cf, while perturbation methods yield a 
maximum near 2kp. The available quantum Monte Carlo data seem to agree with our results. These 
findings imply that many effects, e.g., Friedel oscillations, Kohn anomalies, effective electron-electron 
interactions, etc., are affected by the strong correlations among anti-parallel spins in 2D. 
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Introduction. The study of the uniform two- 
dimensional electron gas (2DEG) continues to bring out 
new and unexpected properties. Its physics depends cru- 
cially on the "coupling parameter" T = (potential en- 
ergy) /(kinetic energy) which defines the strength of the 
Coulomb interactions. The T for the 2DEG at T = and 
mean density n turns out to be equal to the mean-disk 
radius r s = (7m) -1 / 2 per electron. The density parame- 
ter r s , the spin polarization £ and the temperature T are 
the only physical variables that appear in this problem. 
As the coupling constant increases, the 2DEG becomes 
more like a "liquid" , but here we continue to loosely call 
it an electron "gas" . An important property of the 2DEG 
which captures many of the exchange and correlation ef- 
fects and central to its physics is the response function 
x(k,oj). It is expressed in terms of a zeroth-order re- 
sponse x°(fc, uu) and a local-field factor (LFF), also called 
a local-field correction (LFC), denoted by G(k,oj). 

x(k,Lu) = x(k,oj)°/[l - V k (l - G(k,u))) X (k,u)° (1) 

The LFF is also closely related to the vertex function 
A(k,u>). The simplest static form, G(k), is identical with 
G(k, 0) at k = 0, and begins to differ from G(k, ui) as k or 
lo increases. In general, the G(k, to) for u> smaller than the 
plasma frequency uo p behaves like a static quantity. Hence 
the use of a static form G(k) is often adequate. As such, 
considerable effort has been devoted to determining G(k), 
using perturbation theory, kinetic equation methods fj\, 

etc. A partially analytic, partially empirical, approach 
is often implemented by invoking parametrized models 
constrained to satisfy sum rules ||. These sum rules 
invoke physical quantities which are usually obtainable 
only from numerical simulations. While these efforts have 
considerably extended the theory beyond the random- 
phase approximation (RPA) response, its validity is still 
restricted to the low-r s regime. 

We have shown in a number of recent papers [S ||, |(| 
that the static properties of the 2D and 3D electron 
systems can be calculated via an equivalent classical 
Coulomb fluid having a temperature T q such that it has 



the same correlation energy as the quantum system at 
the physical temperature T — 0. The "quantum temper- 
ature" T q was shown to be given by 



t = Tg/Ep = 2/[l + 0.86413(r]/ 6 - l) 2 ] 



(2) 



where Ep — 1/r 2 is the Fermi energy in Hartree atomic 
units. At finite temperatures, the classical-fluid temper- 
ature T c f is taken to be (T 2 + T 2 ) 1 / 2 , as discussed in 
Ref. [|. The pair-distribution functions (PDFs) of the 
classical fluid are calculated using the hyper-netted-chain 
(HNC) equation j?J inclusive of bridge corrections. A set 
of three coupled equations, for gn, g i2 and 322 for the 
two spin-species i=l,2 is solved. This classical mapping 
of quantum fluids within the HNC was named the CHNC. 
The static response function and the static structure fac- 
tor of a classical fluid are identical (except for numerical 
factors), and we also showed ||, || how this can be ex- 
ploited for a very simple determination of the static LFF 
of the quantum electron fluid. In ref. f| we showed that 
various sum rules, e.g, the behaviour of G(k) at large k 
and its relation to the "on-top" value g(0) of the pair- 
distribution function (PDF) are satisfied by this "classi- 
cally" determined G(k), and that the results agreed well 
with, e.g., the 3D-LFF of Utsumi and Ichimaru ||. 

The application of the CHNC method to the 2DEG 
required the inclusion of short-range clustering effects 
("bridge terms") going beyond the usual HNC approx- 
imation. These contributions seem to play a role similar 
to the "back-flow" contributions used in QMC simula- 
tions ||. Such terms were necessary for gi2, while the 
need was less stringent for gn since the Pauli exclusion 
blocked the clustering for identical spins. However, it is 
well known that the compressibility sum rule is violated 
by the HNC approximation, and that a bridge term is in 
principle necessary [|l0|, even for gu. While this short- 
coming is not strongly felt in the calculated distribution 
functions |yj, and even less so in the energy, it needs 
to be addressed if the small-fc limit of the LFF is to be 
correctly recovered. 

The objective of this work is to evaluate the LFF of 
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the 2DEG as a function of r s and show that it does not 
behave in the manner expected from standard perturba- 
tion calculations using the RPA screened Coulomb 
interaction. Such calculations give LFFs with a "hump" 
at 2&f due to the singular nature of x{k = 2kp), while 
we find that the interactions have moved the hump to- 
wards ~ 3k p. A limited set of QMC data for the 2-D 
LFF is available in the literature jl3[ [l4| and seems to 
be in agreement with our findings. If this is indeed the 
case, effects like the Kohn anomaly, Friedel oscillations, 
effective attractive electron-electron interactions! 12 etc., 



would be weaker in 2-D. This emphasizes the need for 
further work in this field. 

The local-field factor. The LFF was already defined 
m Eq. 0. Here we are concerned with the static form 
G(k), defined with respect to a reference "zeroth-order" 
response function. It is customary to use the Lindhard 
function x(k)° T for this purpose. However, a number of 
author [|l5| pointed out that another natural choice is to 
use the "density-functional" non-interacting form x{k)° 
containing the occupation numbers corresponding to the 
interacting density. Richardson and Ashcroft(RA) |1(| 
explicitly constructed a set of LFFs based on this x(ky, 
containing the corrected occupation numbers. The cal- 
culation of RA was for the 3DEG, and identified the sim- 
plest conserving set of diagrams in the perturbation ex- 
pansion for the RPA-screened interaction. This set (Fig. 
2 of RA) is the same as selected earlier by Geldart and 
Taylor |17| who calculated the static 3D case. In all cases, 
dynamic or static, the LFF calculated from the diagrams 
is further restructured using quantum Monte-Carlo data 
in several ways. In effect, most of the currently proposed 
LFFs, be they for 2D (e.g., @) or 3D (e.g, §), are in 
some sense only partially analytic since they use the nu- 
merical simulation (QMC) data and fit the behaviour at 
k=Q, k — ► oo and possibly S(k) etc., to form the LFF. 
This truly emphasizes the difficulty and delicateness in- 
volved in the determination of the LFF. 

The CHNC also uses the QMC correlation-energy 
E c (r s ), at the initial stage of constructing the "quantum 
temperature" T q (r s ) of the classical fluid at zero temper- 
ature (T — 0). From then on we have a self-contained 
method for obtaining all the static properties including 
E c (r s X,T) and its gradients, as well as g(r), S(k) etc. 
The CHNC also provides a very simple formula for the 
LFF, via the classical-fluid LFF. Unlike in the quantum 
case, for a classical fluid, x(fc) is directly related to the 
structure factor. 

S ij (k) = -(l/(3) Xij (k)/(n i n j ) 1 / 2 (3) 

Hence, taking the paramagnetic case for simplicity, 



V c (k)G(k) = V c {k) 



n 



S(k) S°{k) 



(4) 



Here T c f equals T q if the physical temperature T = 0. 
In CHNC, and hence in these expressions, the x°(k) and 
5°(fc) are based on a Slater determinant, while the Lind- 



hard function is applicable to the non-interacting case 
without antisymmetrization of the wavefunction. 

The above expressions show that the LFF is immedi- 
ately available if the interacting and noninteracting struc- 
ture factors are known. These are explicitly known since 
the HNC+Bridge equations for the classical fluid yield 
the PDFs gij (r) and Fourier transforms of which directly 
yield the Sij (k) needed here. Alternatively, any other 
source of S(k), e.g., QMC, may be used, while S°(k) for 
the 2DEG is analytically known. Equation |i| involves a 
difference between the inverses of S(k) and 5°(fc), and 
hence it is very desirable to calculate them to the same 
accuracy, since their difference has to cancel the Coulomb 
potential which becomes large as k — > 0. The Coulomb 
potential V c {r) for point-charge electrons is 1/r. How- 
ever, the classical electron at the temperature T c f is lo- 
calized to within a thermal wavelength. Thus, we use the 
"diffraction corrected" form j(| 

V c (r) = {l/r)[l-e- rk ^\ (5) 

V c (k) = - (k\ h + fc 2 )- 1/2 ] (6) 

The k t h was determined by numerically solving the 
Schrodinger equation for a pair of 2-D electrons in the 
potential 1/r and calculating the electron density in each 
normalized state [fl§l ■ It was found that 



hh/k 



°„ = 1.158T c y 03 



where T c f is in a.u., Here k® h is the de Broglie thermal 
wavevector (2irm*T c f) 1 ' 2 , with m* , the reduced mass of 
the scattering electron-pair, equal to 1/2. 

The introduction of the de Broglie wavevector kth adds 
a new fc-scale which competes with the Fermi wavevector 
kF of the non-interacting problem. However, k t h becomes 
equal to about kp only when we approach the Wigner 
crystallization ( r s ~ 35 ) regime. Hence, as far as this 
study is concerned, the critical wavevector k c which sep- 
arates the small-fc region from the large-fc region will be 
taken to be 2fcp . Although Eq. || is already very simple 
to compute, the explicit cancellation of V c by the terms 
in the 1/5 — l/5o can be realized as follows. For the 
paramagnetic case, we have: 

S{k) = I + nxih^k) + h 12 (k)] (7) 
= - nx{ Cll (k) + c 12 (k)}} (8) 

Here we have used the Orstein-Zernike relations and 
the direct (cjj) and total (fojj) correlation functions, 
while a;=l/2 is the fractional composition. The pair- 
distribution function Q in CHNC has the form 

gij(r) — e l~^ i j( r )~ Ci j( r )+' li j( r )+- Bi j( r )] (9) 



<kj(r) = P l3 {r)S tJ +V c (r) 



(10) 



Here Pij is an effective potential for the Pauli-exclusion 
effect, and is such that, when used in the HNC equations, 
reproduces the non- interacting PDF, g°(r). On defining 
modified (short-ranged) direct-correlation functions: 



Cij(r) = Cij(r)+P<l>ij(r) 



(11) 
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we can express the LFF in a form where V c (k) has been 
explicitly cancelled out. 

G(k)=Tx[c n +c 12 +B u (k)+B 12 (k)-c n }/V c (k). (12) 

This is free of numerical inaccuracies at small-k, and con- 
tains the bridge terms that may be used to recover the 
compressibility sum rule etc. On the other hand, this re- 
quires explicit forms for the bridge functions (which can 
be constructed using the hard-disk model bridge func- 
tion ||). However, here we propose a much simpler, 
practical, approach based on using Eq. |i| which requires 
only the S(k) and S°(k) which may be available from 
other sources beside the CHNC. 

The 2-D compressibility sum rule Q leads to the fol- 
lowing small-fc form for the LFF. 



G(k) = C k k 



Co 



1 

7T 



4 





rdEc 

dr s 



d 2 E c 



(13) 
(14) 



Here E c is the correlation energy per particle in atomic 
units, and a = 1/^/2. Several fit formulae (based on 
QMC data, |§ glj) are available for E c (r s ). Hence 
Co can be calculated easily. Equation [l3| is valid only 
for small-fc where G(k) is linear in fc. Based on more 
extensive calculations using Eq. [l^, we assume a linear 
interpolation for the small-fc region fci to fc c , where fci is 
the smallest value of k in our S{k) tabulation. That is, 
we replace the form: 



G{k) = 1 - 



T, 



of 



1 



1 



nV c (k) lS(k) 
by the constrained form: 



G(fci) 
G(k) 



C fci 

G(k c )-G(h) 



G(fc) = 1- 



T, 



1 



S°(k) 



k — fci 
fc c fci 
1 



all k 



(15) 



nV c (k) IS(k) S°(k) 



k < fc c . (16a) 
fc > fc c .(16b) 



We have also displayed the 3DEG local-field factor for 
r s = 1 and 5. The 3D-LFF shows a weak maximum 
near 2fc^ and tends to 1 — g(0) for large fc. The QMC 
data (triangles) are extracted from Fig. 1 of ref. (l|, and 
are not available for r s —2Q. Note that these do not 
tend to the 1 — g(0) but has a lincar-g dependence for 
asymptotically large -q, consistent with a zeroth-order 
Lindhard response. For r s = 1, the thermal wavevec- 
tor (kth/kF = 2.88) is larger than 2kp. Hence the re- 
gion where G(fc) is linear in fc may extend well beyond 
2kp . However, we have retained a linearization up to 
2kp which joins with the fc > 2fc^ behaviour given 
from CHNC. Also, the extended linearity in the QMC 
is partly due to its large-g asymptotic behaviour. For 
r s =5 strong-coupling effects show up and a broad hump- 
like structure near k/kp ~ 3 is manifest. When we 
go to r s = 10 this behaviour is more pronounced, and 
sharpens even more for r s =20. The required large-fc be- 
haviour, G(k) — > 1 — <?(0), of the LFF is automatically 
satisfied by CHNC and hence we do not impose it. The 
difference between 1/S(k) and 1/S°(k) required in Eq. [l5| 
increases for strong coupling, and the direct method of 
calculating the LFF, using Eq. ||, instead of Eq. [ij] be- 
comes adequate. Thus we compare the LFF calculated 
without imposing the compressibility sum rule, and using 
S(k) from two different sources and different (but equiv- 
alent) formulations. For the cases r s = 10 and 20, we 
have given (dashed line) LFFs evaluated purely using the 
S(k) and S°(k) derived from CHNC, and using Eq. [j^. 
These include only the Byi bridge term. The linear form 
based on the compressibility sum rule has not been im- 
posed. We also present the LFFs (solid lines with circles) 
obtained using the S(k) given in 1989 by Tanatar and 
Ceperely|pj||, together with the linearized forms (solid 
lines with squares) resulting from them. These results 
show that the maximum around 3k p is obtained irre- 
spective of the method used, while the small-fc region is 
clearly sensitive to the numerical procedures used. More 
details about the CHNC method and on-line access to 
our codes may be obtained at our website ]22|. 



Here we choose fc c = 2kp. In practice, fc c is chosen 
to be the k-point closest to 2fc^ in the S(k) tabulation 
which gives a continuous accord between the low-fc and 
high-fc regimes, and hence k c /kp may be, e.g., 1.9 or 2.1. 

Results. We now consider the LFF for r s = 1,5,10 
and 20 calculated using Eq. |[ and displayed in Fig. 0. 



conclusion Unusual features of the 2D-LFF not found 
in the 3D-case, and unexpected from perturbation theory, 
have been revealed via the CHNC method. These results 
seem to be in agreement with the limited QMC results 
currently available. They should stimulate further work 
on the important problem of the 2D response using non- 
perturbative methods. 
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